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Abstract. We present new results from the VLBA 2 cm Survey, an imaging survey of active galactic nuclei (AGN) at sub- 
milliarcsecond resolution. We extend the structural variability monitoring program of more than 130 extragalactic parsec-scale 
radio jets from over 170 AGN to a total of eight years. The sample is explored further in time for all sources, and intensively for 
individual objects. We report new detailed results on the compact sources IIIZw2, AO0235+16, and NRA0512. 



1. Introduction 

Over the last decade we have used the Very Long Baseline 
Array (VLBA) at 15 GHz to image the structure of more than 
170 active galactic nuclei (AGN). This large collaboration is 
known as the "2 cm Survey". Most of the sources show a core- 
jet morphology with typical jet structures (e.g., Zensus 1997 1. 
The 2 cm Survey aims to study the jet kinematics in AGN to 
better understand the processes involved in the formation and 
evolution of powerful jets. The first images of 132 sources, 
selected from data taken between August 1994 and October 
1996, were presented by Kellermann et al. ([1998 1. An addi- 



tional 39 sources were presented by Zensus et al. ( 2002 1, from 
data collected from August 1997 to March 2001. Many of the 
jets may be described by a small number of apparently discrete 
features, which are generally called components. Kellermann et 



al. (2004 1 reported the sky motions over the above mentioned 
time lapse of 208 components in 110 AGN from the sample. 
We have made all the survey images available on the Internet 
under http : //www . nrao . edu/2cmsurvey. 

We continued adding new observations from October 2001 
to October 2002. We also included data from other survey 
projects in our kinematic analysis, such as the pre-VSOP ob- 
servations from L. I. Gurvits et al. (in prep.), the S5 Polar 
Cap Sample observations by Perez-Torres et al. (2004i, and 
the observations from GHz-peaked-spectrum radio sources 
by Stanghellini et al. (2001 1. To continue beyond the obser- 
vations reported here, we have assembled a complete flux- 
density limited sub-sample from the full 2 cm VLBA sur- 
vey. The selection criteria are described in Kellermann et 



al. (2004 1. Observations of the complete sample are being 
continued since May 2001, under the name of " Mo nitoring 
of Jets in AGN with VLBA Experiments (MOJAVE)", see 
http : //www . physics . purdue . edu/astro/MO J AVE/. This 
program includes also dual polarisation observations. 

In this contribution we present results based on the 2001- 
2002 observations, which turn out to be complementary to the 
previously published ones. 



2. Kinematics 

The survey observations reported so far (Kellermann et al. 
1998 Zensus et al. |2"0"0"21 Kellermann et al. 2004 1 included 



29 observing sessions. We added 18 observing sessions with 
typically eight sources per session, where each source was ob- 
served once per hour for 6 min to 8 min over a range of 8 h in 
hour angle. The data reduction procedures for imaging, com- 
ponent identification and measurement of its position, combin- 
ing difmap and JilVS tasks, are described in Kellermann et 
al. (Q998; 2004 1. Our additional data allowed us to analyse the 
kinematics of 41 components in 19 sources. The measured po- 
sitions and speeds are presented in Fig.[2and tabulated column 
5 of Table ^ in some cases, when the jet features were not 
well resolved, or in other peculiar cases, we performed a fit to 
the visibilities using Gaussian functions. In this case, we have 
marked an asterisk at the name of the sources in Fig.^ 

Apart from the 110 sources reported in Kellermann et al. 
(2004) and the 19 ones in this paper, some objects remain 
without kinematic analysis. Those objects are listed in Tabled 
From the remaining objects, some of them have such a structure 
that prevents a reliable analysis (e.g., 0521-365 hosts a long jet 
and the low declination makes mapping and model fitting to be 
difficult). A fraction of the objects have a parsec-scale jet where 
kinematics can be measured, but only two observing epochs 
have been observed so far. Some of those 1 are being monitored 
at the MOJAVE program, and therefore reliable kinematics will 
be available in the near future. The rest 2 are not being observed 
further, so unless alternative observations become available, no 
kinematic study using more than the existing two epochs will 
be possible. Finally, the objects not categorised above are com- 
pact, and some detailed (u, v)-model fitting can eventually re- 
veal changes in the structure, but higher frequency observations 
should provide the final answer on their kinematics. 

The complete set of kinematics for the 110 sources reported 
in Kellermann et al. (2004i and the newer ones, together with 
an analysis of the component (x, y) motions and the flux density 
evolution in selected cases will be reported in Ros et al. (in 
prep.). 



1 0119+115, 1504-167, 1741-038 (OT-068), and 2155-152 
(OX -192). 

2 0108+388 (OC+314), 0201+113 (OD+101), 0405-385, 
0723-008 (OJ+039), 0859+470 (4C +74.29), 1511-100 (OR -118), 
1514-241 (AP Librae), and 2255-282. 
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Fig. 1. Plots showing the change in separation with time of features of one-sided jets for which we have measured a velocity from 
observations at three or more epochs, after including the data from the newer observations and from another 15 GHz surveys (see 



text). The codes used are the same as in Fig. 1 of Kellermann et al. (2004i. An asterisk denotes sources that were model fitted 
in the (u, v)-plane rather than in the image plane. The solid lines denote the best least-squares linear fit to the data, and the slope 
represents the proper motion, fi. 



3. Compact, but slightly resolved sources 

We have reported in detail in Sect. |2] the kinematics of the 
sources with a core-jet morphology. In our sample we also find 
many compact objects. A careful look to these objects may 
yield surprises, since not all of them are unresolved. Here we 
show three typical cases of this kind. An analysis of the fine 
scale structure based on the interferometric visibilities is pre- 
sented in Y. Y. Kovalev et al. (in prep.). 

IIIZw2 (0007+106)— cracking the nut ? The first panel 
in Fig. H reports the proper motions IIIZw2 derived from a 
model fitting of the interferometric visibilities using two circu- 
lar Gaussian components — case (i); the use of two delta dis- 
tributions (case (ii)) yields equivalent results. This Seyfert 1 
galaxy, also classified as "radio intermediate quasar", is highly 
variable, and a radio flare (Falcke et al. 1999 1 happened in the 
course of our VLBI monitoring, at epoch ~1998.8. Brunthaler 
et al. (2003 1 analysed the structure of the radio source at 
15 GHz model fitting the visibilities with two point components 
in nine epochs from 1998.2 to 2000.7. Their data show a con- 
stant change in the separation of the components with speeds of 



~ 0.6 c. A direct comparison of their results and our plot show 
discrepancies, since our model at epoch 1998.21 shows a larger 
component separation, although our 1998.83 and 2002.02 data 
points are compatible with their measurements. A further ap- 
proach to understand the changes in the source is to model fit 
its structure with a single circular (case (iii)) or elliptical (axis 
ratio 0.4, case (iv)) Gaussian component. The flux densities and 
the corresponding brightness temperatures are shown in Fig.|2] 
Notice that the brightness temperature was stable while the flux 
density was rising, and remained low at the highest state of the 
flare, becoming much higher in the post-flare phase and coming 
slowly back to the values of the rising phase at the last epoch. 
We can compare this result with the suggestions of Brunthaler 
et al. (2003 1 that the source can be explained with an "expand- 
ing balloon" model, say, "a radio galaxy in a nutshell". Our 
comparatively low brightness density for the rising and peaking 
phase of the flare is in apparent contradiction with the interpre- 
tation of the 15 GHz data being the post-shock material from 
the hot-spots expanding in ultra-compact scales and interacting 
with the environment. In any case, we observe the same pro- 
cess of inflation and deflation in the source, but with a sparser 
time sampling and at only one frequency. 
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Table 2. Sources without reported kinematics 



Table 1. Source Kinematics 



Object 


Optical 

class 3 


Id 


N 


H 

[mas yr -1 1 




to 
[yd 


0007+106 


G 


B 


6 


-0.011 ±0.015 


-0.06 ± 0.09 




0221+067 


Q 


B 


3 


0.173 +0.013 


5.26 ± 0.41 


1963.0 


0834-201 


Q 


B 


3 


-0.021 ± 0.005 


-2.0 ±0.5 




0836+710 d 


Q 


B 


3 


-0.03 ± 0.08 


-2.6 ±6.7 








C 


3 


0.118 ±0.017 


10.1 ± 1.4 


1994.2 






D 


3 


0.105 ±0.021 


9.0 ± 1.8 


1985.1 






E 


3 


0.096 ± 0.022 


8.3 ± 1.9 


1970.0 


1145-071 


Q 


B 


3 


0.039 ± 0.005 


2.52 ± 0.34 




1148-001 


Q 


B 


3 


0.089 ± 0.004 


7.23 ± 0.32 


1978.5 






C 


3 


0.115 ±0.003 


9.32 ± 0.24 


1958.8 






D 


3 


0.031 ± 0.007 


2.53 ± 0.60 








E 


3 


0.118 ±0.029 


9.7 ± 2.3 




1354+196 6 


Q 


B 


3 


0.115 ± 0.028 


4.7 ± 1.2 


1990.2 






D 


3 


0.202 ± 0.025 


8.2 ± 1.0 


1973.5 


1458+718 6 


Q 


B 


5 


0.102 ±0.012 


4.9 ± 0.6 


1993.8 






C 


5 


0.11 ±0.06 


5.3 ± 2.9 








D 


5 


0.010 ± 0.033 


0.5 ± 1.6 




1502+106 


Q 


B 


3 


-0.02 ± 0.06 


-1.4 ± 4.3 








C 


3 


0.11 ±0.05 


8.7 ±4.2 


1985.1 


1504+377° 


G 


B 


3 


0.032 ± 0.023 


1.2 ±0.9 


1989.6 






C 


3 


0.107 ± 0.006 


4.11 ± 0.25 


1992.2 






D 


3 


0.156 ±0.008 


5.98 ± 0.30 


1990.5 






E 


3 


0.19 ±0.09 


7.3 ± 3.4 


1981.1 


1538+149 


Q 


B 


4 


-0.003 ± 0.003 


-0.11 ±0.10 








C 


4 


-0.033 ±0.021 


-1.2 ±0.7 








D 


4 


0.027 ± 0.027 


0.97 ± 0.94 




1555+001 


Q 


C 


3 


0.08 ± 0.07 


6.0 ±5.6 




1622-253 


Q 


B 


3 


-0.059 ± 0.022 


-2.5 ±0.9 




1638+398 


Q 


B 


5 


0.003 ± 0.008 


0.24 ± 0.57 








C 


5 


0.139 ±0.008 


10.2 ±0.6 




1954+513 


Q 


B 


3 


0.22 ± 0.09 


13.2 ±5.1 


1993.2 






c 


3 


0.28 ± 0.13 


17.1 ± 7.9 








D 


3 


-0.03 ± 0.07 


-1.6 ±4.3 








E 


3 


0.065 ± 0.009 


3.9 + 0.6 








F 


3 


-0.08 ± 0.07 


-4.5 ±4.5 








G 


3 


-0.036 ± 0.019 


-2.19 ± 1.17 




2121+053 


Q 


B 


3 


0.071 ±0.006 


5.64 ± 0.46 


1993.2 


2128+048 f 


G 


B 


3 


-0.007 ± 0.030 


-0.37 ± 1.55 








C 


3 


0.003 ± 0.001 


0.16 ±0.04 




2227-088 c 


Q 


B 


3 


0.02 ±0.16 


1.3 ± 11.1 




2318+049 6 


Q 


B 


5 


0.170 ±0.042 


6.1 ± 1.5 


1970.7 



a Optical classification according to the Veron-Cetty & Veron 2001 i catalog, where 

Q = quasar and G = galaxy. 
b To compute the linear speed, and throughout the paper we use a cosmology with 

H = TOkms- 1 Mpc -1 ,^ = 0.3, andH A = 0.7. 
Extrapolated epoch of origin. 

d Includes data from the S5 Polar Cap Sample (Perez-Torres et al. 2004 >. 
L " Includes data from the pre-VSOP observations by L. I. Gurvits et al. (in prep.). 
f Includes data from observations taken in April 1996 on Gigahertz-peaked sources and 
published by Stanghellini et al. |2001| . 



AO 0235+1 6 (0235+164) & NRA0512 (1638+388) Those 
two AGN show compact structure (left panel in Fig. [3} and 
flux density variability. The BLLac object AO 0235+16 shows 
intra-day variability (Quirrenbach et al. 1992 Kraus et al. 
1999 1 and recurrent outbursts with a period of ~(5.7 ± 0.5) yr 
(Raitieri et al. 2001). It is well known that a very high bright- 
ness temperature component is present at milliarcsecond-scales 



Optical 



UDject 


class 3 


IN 


Comment 


MOJAVE b 


0048-097 


B 


5 


Compact, N emission 


Y 


0108+388 


G 


2 


EW structure 


N 


0119+115 


Q 


2 


N jet 


Y 


0201+113 


Q 


2 


NWjet 


N 


0218+357 


Q 


4 


Gravitational lens c 


N 


0235+164 


B 


5 


Compact, E emission 


Y 


0310+013 


Q 


1 


Compact, SE emission 


N 


0405-385 d 


Q 


2 


Wjet 


N 


0420+022 


Q 


2 


Compact 


N 


0521-365 d 


G 


3 


Long NW jet, difficult comp. id. 


N 


0552+398 d 


Q 


4 


Compact, W emission 


Y 


0602+673 


Q 


4 


Compact 


N 


0615+820 d ' e 


Q 


7 


Compact, halo 


N 


0723-008 


B 


2 


NWjet 


N 


0808+019 


B 


4 


Compact, S emission 


Y 


0859+470 


Q 


2 


Njet 


N 


1032-199 


Q 


2 


SW emission 


N 


1124-186 


Q 


2 


Compact, S emission 


Y 


1155+251 


Q 


4 


Too complex 


N 


1328+254 


Q 


4 


Too complex 


N 


1328+307 


Q 


4 


Too complex 


N 


1354-152 


Q 


2 


Compact, NE emission 


N 


1424+366 d 


Q 


5 


Compact, SW emission 


N 


1504-167 


Q 


2 


SEjet 


Y 


1511-100 


Q 


2 


Ejet 


N 


1514-241 


B 


2 


S jet 


N 


1519-273 


B 


2 


Compact, W emission 


N 


1739+522 


Q 


2 


Compact, NE emission 


Y 


1741-038 


Q 


2 


SW emission 


Y 


2155-152 


Q 


2 


SWjet 


Y 


2255-282 


Q 


2 


SWjet 


N 



3 Optical classification according to the Veron-Cetty & Veron > 2001 > catalog, where 

Q = quasar, B = BL Lac object, and G = galaxy. 
b Belonging to the complete sub-sample (MOJAVE, observed from May 2002 on). 
e A detailed study of the relative motions within the sub-images A and B and between 

them will be presented elsewhere. 
d Includes data from the pre-VSOP observations by L. I. Gurvits et al. (in prep.). 
e Includes data from the S5 Polar Cap Sample (Perez-Tones et al . 120041 . 



(e.g., Frey et al. EUHOl Peng & de Bruyn 2004 1. The quasar 
NRA0512 shows unresolved structure at 8.4 GHz or 2.3 GHz 
(e.g., Fey & Chariot 2000). For this reason, it has been tra- 
ditionally used as a phase-reference calibrator for the neigh- 
bour source 3C 345 (e.g., Bartel et al. 1986). The clean images 
shown at the left panels of Fig. [3] shows that both NRA0 512 
and AO 0235+16 are compact objects but a detailed analysis 
shows that those are not unresolved. Figure^shows measured 
kinematics for NRAO 512 ussing a three-Gaussian model. The 
multi-epoch model fitting with three components appears con- 
sistent. 

A rough analysis, presented in Figure [3] shows that a sin- 
gle Gaussian fit to the structure of those two objects does not 
adequately reproduce their structure. The mid panels show the 
residual maps after this fitting. The right panels show the resid- 
uals of this fit, and the discrepancy beyond 300 MA is obvious, 
especially in the case of NRAO 5 12. A model with a clean com- 



Eduardo Ros: Extending and Exploring the 2 cm Survey Sample 



III Zw 2 (0007+106) 




1996 



1998 
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Peng, B., & de Bruyn, A. G. 2004, ApJ, 610, 151 

Perez-Torres, M. A., Marcaide, J. M., Guirado, J. C, & Ros, E. 2004, 

A&A, in press |arXiv:astro-ph/0408581 1 
Quirrenbach, A., Witzel, A., Krichbaum, T. P., et al. 1992, A&A, 258, 
Raiteri, C. M„ Villata, M., Aller, H. D., et al. 2001, A&A, 377, 396 
Stanghellini, C, Dallacasa, D., O'Dea, C. P., et al. 2001, A&A, 377, 

377 

Veron-Cetty, M.-P, & Veron P. 2001, A&A, 374, 92 
Zensus, J. A. 1997, ARA&A, 35, 607 

Zensus, J. A., Ros, E., Kellermann, K. I., et al. 2002, AJ, 124, 662 



Epoch (yr) 

Fig. 2. Time evolution of the total flux density and brightness 
temperature in III Zw 2, derived from our model fitting (cases 
(iii) and (iv)). The error bars are formal standard errors. 



ponent and an extended structure can better fit the data — a sort 
of "prussian hat". An analysis of the change of the flux den- 
sities of both the ultra-compact and the extended component 
over time will be presented elsewhere. As in other examples, 
the ultra-compact component is needed to explain the intra-day 
variability (either if intrinsic or if caused by inter-stellar scat- 
tering), and the extended emission is related to a jet structure, 
that can be more or less weak, and face-on or slightly tilted and 
therefore present in the images. 

The analysis of those sources and the rest of the 2 cm 
Survey will keep providing interesting results and shedding 
more light on the nature of the AGN in radio frequencies and 
its relationship with their properties in other wavebands (see, 
e.g., Kadler et al., these proceedings, concerning X-ray obser- 
vations). 
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Fig. 3. Images and visibility residuals from the compact, but resolved sources AO 0235 + 1 6 and NRAO 5 1 2 at their first observing 
epochs in the 2 cm Survey after fitting them with a single Gaussian component. Left panels: hybrid maps. Central panels: residual 
maps after fitting a single elliptical Gaussian component to the visibilities (with flux [mJy], major axis [mas], axis ratio, and 
position angle [deg], respectively of 766:0. 357:0. 55:-l 1.4 and 914:0.271:0.924:24.0). The profile of the ellipse is shown at the 
centre, being the two major axes the full-width half maxima. Notice the need for extended emission beyond one beam size, shown 
in gray scale. Right panels: residuals of the visibility amplitude in Jy versus the interferometric distance, showing the result of 
subtracting the Gaussian model fit from the observed data. 



